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SECTION  I 


INTRODUCTION 

In  this  report  the  integral  equation  for  the  linearized 
supersonic  unsteady  potential  flow  over  a  wing  will  be  derived 
and  the  special  properties  of  its  solutions  will  be  discussed  by 
means  of  limiting  cases  that  can  be  treated  analytically.  As  in 
every  formulation,  by  means  of  an  integral  equation,  the  basic 
prerequisite  is  a  fundamental  solution.  The  expression  from 
which  we  start  here  represents  a  local  source  with  arbitrary  time 
dependence  in  a  system  of  coordinates  fixed  with  respect  to  the 
wing.  (The  undisturbed  air  moves  with  respect  to  this  system  of 
coordinates  from  left  to  right.)  Only  solutions  representing 
outgoing  waves  are  admitted.  The  fundamental  solution  used  in 
the  integral  equation  represents  a  time  dependent  doublet;  it  is 
obtained  from  the  expression  for  a  source  by  a  differentiation 
into  the  direction  normal  to  the  planform.  Use  of  such  a 
formulation  in  subsonic  flows  is  shown  in  References  1  and  2. 

The  corresponding  solutions  for  supersonic  flow  are  given  by  the 
same  formula.  Nevertheless,  an  ambiguity  aries.  It  seems  as  if 
one  has  two  instead  of  one  fundamental  solution,  each  of  which 
satisfies  the  flow  differential  equation  in  the  region  in  which 
the  derivatives  exist.  At  this  stage,  the  typical  phenomena  of 
supersonic  flow  are  encountered.  One  has  regions  of  influence, 
with  a  boundary  given  by  the  surface  of  the  Mach  cone  through  the 
point  where  the  source  is  located.  Along  this  surface  the 
fundamental  solutions  are  singular.  Outside  the  Mach  cone  they 
are  identically  equal  to  zero.  One  might  try  to  resolve  the 
above  mentioned  ambiguity  by  discussing  whether  the  partial 
differential  equation  for  the  potential  is  satisfied  in  the  sense 
of  generalized  functions  as  one  passes  through  the  surface  of  the 
Mach  cone.  Some  ideas  which  lie  in  this  direction  are  found  in 
Appendix  C.  The  ambiguity  is  resolved  in  Section  II. 

The  point  of  departure  is  the  fundamental  solution  for  a 
time  dependent  local  source  in  a  system  of  coordinates  at  rest 
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with  respect  to  the  unperturbed  air.  In  this  formulation,  no 
ambiguity  is  encountered.  Next  one  proceeds,  still  in  the  same 
coordinate  system,  to  a  distribution  of  such  sources  along  the  x- 
axis;  and  subsequently,  by  a  specialization  to  a  localized  time 
dependent  source  moving  along  the  x-axis  from  right  to  left  with 
the  desired  free-stream  velocity.  Finally,  a  transformation  to 


coordinates  in  which  the  undisturbed  air  moves  with  the  desired 
free  stream  velocity  from  left  to  right  is  carried  out.  This 
gives  the  desired  fundamental  solution  (a  local  time  dependent 
source  in  this  last  system  of  coordinates).  This  procedure  is 
applicable  for  subsonic  as  well  as  supersonic  flows.  In  subsonic 
flow  one  finds,  of  course,  the  particular  solutions  used 
previously.  In  supersonic  flow  one  finds  the  sum  of  two 
expressions,  each  of  which  satisfies  the  flow  differential 
equation  where  the  derivatives  exist.  The  derivation  gives  only 
one  linear  combination  of  these  expressions  and  thus  removes  the 
ambiguity . 

In  Section  III,  familiar  formulae  expressing  the  upwash  in 


terms  of  the  motion  and/or  deformations  of  the  wing  are  derived. 

Section  IV  formulates  the  integral  equation.  The  flow  field 
is  represented  by  time  dependent  doublets  distributed  over  the 
planform.  The  source  strength  in  its  dependence  upon  time  and 
location  on  the  planform  is  unknown  at  this  stage.  An  essential 
part  in  the  formulation  of  the  integral  equation  is  a  limiting 
process  by  which  one  determines  the  upwash  due  to  this  doublet 
distribution  as  one  approaches  the  planform  from  above  or  below. 


(The  results  are  the  same.)  In  the  final  formulation,  this 
limiting  process  no  longer  appears. 

In  the  actual  computation,  the  upwash  due  to  the  motions 
and/or  deformations  of  the  wing  is  built  up  from  expressions  each 
of  which  is  the  product  of  some  function  of  coordinates  within 
the  planform  and  a  Hamilton  step  function  in  time.  It  is  then 
sufficient,  if  one  solves  the  integral  equation  (separately  for 
each  function  of  the  planform  coordinates  for  this  special  time 
dependence).  For  such  an  upwash,  the  flow  tends  towards  a  steady 


state.  In  contrast  to  subsonic  flows,  the  steady  state  is 
reached  within  a  finite  time.  The  essential  physical  phenomena 
occur  during  this  transition  time. 

Section  V  deals  with  two  limiting  cases  that  can  be  treated 
analytically.  The  first  one,  the  two-dimensional  steady  flow,  is 
very  simple;  the  contributions  of  the  integrals  in  the  integral 
equation  are  zero.  The  second  example  is  the  one-dimensional 
unsteady  flow.  Again,  the  given  upwash  is  built  up  by  step 
functions  in  time.  (There  is  no  dependence  upon  the  point  of  the 
planform  under  consideration.)  For  a  step  function,  one  can 
determine  the  region  in  which  the  integrand  is  different  from 
zero.  Within  these  regions  the  iterations  can  be  carried  out  in 
a  closed  form,  and  one  obtains  the  expected  result.  For  a 
general  upwash  (but  still  with  the  time  dependence  given  by  a 
step  function),  the  regions  of  integration  are  the  same,  except 
that  there  are  no  doublets  upstream  of  the  leading  edge.  The 
example  shows  that  the  integrands  jump  from  a  finite  value  to 
zero,  as  one  passes  its  integrands  over  the  boundary  of  the 
region  integration.  The  boundaries  of  the  regions  of  integration 
do  not  coincide  with  the  subdivisions  of  the  planforms  imposed  by 
a  discretization.  In  an  actual  program,  one  must  therefore 
decide  whether  to  allow  these  jumps  to  be  smoothed  out,  or 
whether  it  is  preferable  to  identify  these  boundaries. 

In  the  initial  stage  of  numerical  experimentation  with  a 
program,  examples  in  which  analytic  solutions  are  available  are 
useful.  In  Section  VI  an  exact  solution  for  the  two-dimensional 
problem  with  a  straight  leading  edge  normal  to  the  flow  direction 
is  derived.  This  is  possible  because  the  flow  field  can  then  be 
represented  by  a  source  distribution  (rather  than  by  a  doublet 
distribution),  and  because  the  strength  of  the  sources  is 
directly  given  by  the  prescribed  upwash.  Again,  a  step  function 
gives  the  time  dependence  of  the  upwash.  If  the  distribution  of 
the  upwash  over  the  planform  is  sufficiently  simple,  for 
instance  given  by  polynomials  in  the  coordinates,  then  the 
integrals  that  arise  can  be  evaluated  analytically.  In  testing  a 


numerical  procedure,  one  will  solve  the  same  problem  by  means  of 
the  integral  equations  for  the  doublet  distribution  and  examine 
how  well  one  recovers  the  analytical  results  for  the  potential. 

In  Appendix  A  certain  integrals  are  evaluated  analytically. 
Appendix  B  determines  the  shape  of  the  regions  of  integration,  if 
one  evaluates  the  integrals  at  a  distance  from  the  planform. 
Originally,  this  was  done  for  the  purpose  of  an  extended 
discussion  of  the  one -d i mens ional  unsteady  problem.  (The 
treatment  in  Section  V  is  restricted  to  points  of  the  planforms.) 
The  shape  of  the  regions  of  integration  obtained  in  this  Appendix 
can  hardly  be  foreseen  on  the  basis  of  the  results  of  Section  V. 
Again,  the  time  dependence  is  given  by  a  step  function.  If  the 
time  after  the  step  is  small,  then  it  is  insufficient  for  the 
perturbations  to  propagate  the  point  off  the  planform  under 
consideration;  regions  of  integration  do  not  exist.  If  the  time 
is  slightly  larger  than  a  certain  cutoff  value,  the  region  of 
integration  is  a  circle  for  one  of  the  expressions  in  the  upwash 
integral;  for  the  other  expression  there  is  no  region  of 
integration.  Regions  of  integration  of  this  type  are  not  present 
if  one  evaluates  the  upwash  at  the  planform.  Beyond  a  further 
cutoff  point,  one  finds  nonvanishing  regions  of  integration  to 
both  terms  in  the  upwash  integral.  Their  boundaries  are  formed 
by  a  nyperbola  and  parts  of  a  circle.  These  regions  deform 
steadily  into  those  which  one  encounters  if  the  upwash  is 
evaluated  at  the  planform. 


SECTION  II 

THE  FUNDAMENTAL  SOLUTION 

The  point  of  departure  of  this  report  is  a  fundamental 
solution  which  represents  a  source  located  at  a  fixed  point  in  a 
system  of  coordinates  with  respect  to  which  the  air  moves  with 
the  constant  velocity  U.  The  strength  of  the  source  depends  upon 
time.  In  the  integral  equation  from  which  the  flow  field  is 
determined,  one  needs  the  expression  for  a  doublet.  It  is 
obtained  from  that  for  a  source  by  a  differentiation.  In  the 
subsonic  case,  the  expression  for  the  source  can  be  derived 
directly  from  the  partial  differential  equation  for  the  potential 
linearized  for  small  deviations  from  a  parallel  flow  with  this 
velocity  U.  It  is  not  entirely  obvious  how  to  make  the 
transition  to  supersonic  flows;  in  fact,  in  an  initial  attempt 

the  author  was  led  to  an  erroneous  result.  A  clear  cut  answer  is 

obtained  if  one  starts  with  a  particular  solution  which 
represents  a  time-dependent  source  at  a  fixed  point  in  a 

coordinate  system  which  is  at  rest  with  respect  to  the  surround¬ 

ing  air,  and  subsequently  derives  from  it  the  expression  in  a 
system  which  moves  with  respect  to  the  air.  Similar  approaches 
can  probably  be  found  in  the  literature.  The  derivation  is  shown 
here  in  order  to  make  this  report  self-contained. 

Let  x,  y,  z  be  a  Cartesian  system  to  coordinates,  t  the 
time,  a  the  velocity  of  sound  in  the  undisturbed  field,  L  a 
characteristic  length,  and  $(x,y,z,t)  the  potential.  The 
differential  equation  for  the  linearized  potential  in  air  at  rest 
is  then  given  by 


Introducing  dimensionless  quantities 

x  =  x/L,  y  -  y/L,  z  -  z/L,  t  -  at/L 


<f>  ( x ,  y ,  z ,  t ) 


a  1  L  1  <j>  ( x  ,  y ,  z  ,  t ) 


where  one  expresses  x,y,z,  and  t  on  the  right-hand  side  by 
x,y,z  and  t,  one  obtains 

0  +4)  ■*■<)>  0  =  0 

xx  vy  zz  tt 

Hence,  for  a  perturbation  with  spherical  symmetry 


rr  r  r 


-  4>  =0 

tt 


where 


,  ~ 2  -2  -2,1/2 
r  =  ( x  +  y  +  z  ) 


The  particular  solution  for  an  outgoing  spherical  wave  is  then 
given  by 

t  4  f(t  -  r)  (4) 

r 

Here  f  is  an  arbitrary  function  of  its  argument.  The  potential 
equation  is  an  expression  for  conservation  of  mass.  No  mass 
sources  arise  for  r  *  0,  for  there  the  partial  differential 
equation  is  satisfied,  but  at  r  =  0,  mass  emerges.  This  is  seen 
in  the  following  manner.  The  expression  (4)  gives  a  radial 
veloc i ty 

-  -  — „f ( t-r )  -  1  f'(t-r) 


where  f'  is  the  derivative  of  f  with  respect  to  its  argument.  Th 
mass  flow  through  a  surface  with  radius  r  is  then 

- 4n p [ f ( t-r )  -  rf'(t-r)] 


where  p  is  the  density.  Then  for  r  =  0,  one  has 


mass  flow  =  -4upf(t). 

Except  for  the  factor  -4irp,  which  will  be  unessential  in  future 
developments,  the  strength  of  the  source  at  time  t  is  given  by 
f  (  t)  . 

Next,  we  consider  a  distribution  of  such  sources  along  the  x 

axis.  One  then  obtains  an  axisymmetric  flow  field.  It  is 

sufficient  to  consider  the  flow  in  a  meridian  plane  with 

coordinates  x  and  y  where  y  is  the  distance  from  the  axis  of 

~  2  ~  2  1  /  2 

symmetry.  Later,  y  will  be  replaced  by  ( y  +  z  ) 

For  a  source  at  the  point  x  =  y  =  0,  one  has 

rU,x,y)  =  [(x-02+  y2]1/2  (5) 

Let  the  time  dependent  source  strength  per  unit  length  of  the  x 
axis  at  the  station  x  =  £  be  given  by  f(t,£)  (an  unessential  factor 
—  4 tt p  is  omitted).  Then, 


♦  U.y.t)  -[  LilzsSLLaldl  al  (6) 

J  rU.x.y) 

This  expression  is  now  specialized  to  a  source  moving  in  the 
negative  x  direction  with  the  velocity  U.  In  the  xt  system,  it 
then  moves  with  the  dimensionless  velocity  M  =  U/a.  As  a  first 
step,  we  choose 


f(t  ,5)  =  gU)6(  t-t) 

where  6(t-i)  i3  the  Direct  delta  function.  Considering  t  as  a 
function  of  one  obtains  a  source  which  wanders  along  the  £ 
axis  and  has  strength  g ( C )  at  time  t  =  t(E)-  No  mass  is  expelled 
at  this  station  for  t  /  x(E).  The  mass  expelled  in  an  interval 
t(0~e  <  t  <  t(C)  +  c  with  c  arbitrarily  small,  is  g(^).  Now  we 
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ft s 


choose  t ( C )  so  that  in  a  coordinate  system,  moving  to  the  left 
with  respect  to  the  x,y  system  with  the  dimensionless  velocity  M , 
the  source  is  always  at  the  origin. 


Then 


5  =  -Mt, 


t  U )  -  -m-1£ 


Substituting  this  into  Eq .  (6),  one  obtains 


<t>  (  x  ,  y ,  t )  = 


f  g( l )6 ( t  +  M  1 £-r( l , x , y ) ) dg 


r(l,x,y) 


To  utilize  the  defining  properties  of  the  delta  function,  we 
choose,  instead  of  £ ,  the  argument  of  the  delta  function  as 
variable  of  integration 


t  +  M~  l  -  r ( E , x , y )  =  v 


[m_1  +  ]d£  =  dv 


d^/dv  = 


M( x-£ ) +  r 


Eq.  (8)  now  assumes  the  form 


4>  (x,y,  t) 


f  Mg(Q6  ( v )  dv 


M( x -£  )  +  r 


here,  \  is  regarded  as  a  function  of  t,x,y,  and  v.  We  determine 
the  values  of  £  for  which  the  argument  of  the  delta  function, 
namely  v,  vanishes.  Then  from  the  first  of  Eqs.  (9) 


It  is  convenient  to  introduce 


l,  *  Mt  -  f, 

x  +  M  t  =  x  (12) 

y  -  y 

One  then  obtains  from  Eq .  (11) 

r  =  £/M  (13) 

Since  r  is  positive,  it  follows  from  the  last  equation  that  only 
such  values  of  i  are  admissible  in  Eq.  (10),  for  which  £  is 
positive.  After  expressing  r  in  terms  of  £,  x,  and  y,  one 
obtains  from  Eq.  (13) 

t,  -  M[  (x-rj£-  ♦  y“  ]  =  0. 

This  quadratic  equation  for  f  has  the  solutions 
o-i  ?  2  2  I/O 

f,  .  =  M(1-M‘)  [-Mx  *  (x  *  (1-M  )y  )  ]  ;  i  =  1,2  ( 1  4 ) 

One  has  for  M  <  1 

?  2  2  1  /  2 

j  L  x  +  (1-M  ) y  ‘  j  “ | > | x |  . 

The  sign  of  the  expression  within  the  brackets  in  Eq .  (14)  is 
therefore  determined  by  the  sign  of  the  square  root.  But,  we  had 
observed  above  that  only  positive  values  of  ?,  give  positive 
values  of  r.  If  M  <  1,  only  the  positive  sign  of  the  square  root 
is  applicable. 


For  M  >  1 ,  one  obtains  real  values  of  f  only  if  the  radicand 
i s  pos itive,  i.e.,  for 

x2-  ( M2 -  1 ) y 2  >  0. 

The  present  discussion  refers  to  the  points  for  which  the 
argument  of  the  delta  function  vanishes,  i.e.,  for  which  the 
contributions  to  the  potential,  in  Eq.  (10),  is  different  from 
zero.  The  last  equation  defines  a  forecone  (x<0)  and  an 
af terconc  (  x>0  5  of  the  point  x-G  ,  v *0 .  As  in  steady  supersonic 
flows  the  potential  of  a  source  at  x  =  0 ,  y  =>:)  is  different  from 
zero  only  within  the  fore-  and  aftercone. 

We  write  the  expression 

-  M(M2-1)~1[Mx  ♦  (X2  -  (M2-1  )y'-)1/?J  (15) 

For  M  >  1,  the  sign  of  the  expression  within  the  braces  is 
determined  by  the  sign  of  Mx .  As  f,  is  positive,  only  positive 
values  of  x  are  admitted.  This  restricts  the  values  of  x  and  y 
to  the  aftercone  of  the  origin  in  the  xy-system.  Both  signs  of 
the  square  root  are  admissible. 

After  introducing  x  and  y  and  expressing  r  by  Eq.  (13),  the 
denominator  of  the  integrand  in  Eq.  (10),  evaluated  for 
v  =  0,  is  found  to  be 

M(x-£)  +  r  =  Mx  -  M£,^  +  5^/M 

*  Mx  +  (1  -M2)M-1£.  i  =  1,2. 

Substituting  5  ^ ,  Eq.  ( 1 H)  ,  one  obtains  (both  for  subsonic  and 
supersonic  flow)  for  v  «  0, 

M  ( x  -  £ )  +  r  =  ±  [x2+  (1-M2)yc]1//t_ 

With  this  result,  one  obtains  furthermore  from  Eq .  (9) 


dv 


r  2 
±[X  + 


, ,  m2,  2.1/2 

( 1 -M  ) y  ] 


v=0 


This  shows  that  at  the  point  where  v  =  0,  (,  will  increase  or 
decrease,  depending  upon  the  sign  of  the  square  root. 

Now  one  can  evaluate  the  integral  in  Eq .  (10)  on  the  basis 
of  the  definitions  of  the  delta  function 

+  £ 

|  <5  (  v )  d  v  =  1 
-e 

<S  (  v )  =  0  ,  v  *  0 


One  obtains 


Mg(^) 

♦<hy.t)  -  I  n -m2),?]'75 

Here,  according  to  Eq .  (12), 


*i  "  *1  "  Mt 

The  values  of  are  obtained  from  Eq .  (14);  only  positive  values 
of  are  admitted.  We  have  found  that  for  M  <  1,  only  the 
positive  sign  of  the  square  root  in  Eq.  (14)  is  admitted;  for 
M  >  1 ,  both  signs  will  occur.  One  thus  has  for  M  <  1 


<J>  (  x  ,y ,  t) 


M  g( -Mt+M( 1 -M2) 1 ( -Mx+ (x2+( 1 -M2)y2) 1 72 ) ] 
r  2  ,,  2 .  2.1/2 

[x  ♦  (  1  -M  )y  J 


(16a) 


and  for  M  >  1 


<P(x,y,  t) 


M 


7~2  777  .7  2.1/2 

[x  -  (M  -1 )y  J 

{  g (  -Mt  -  M(M2-1)  1 [ -Mx+ (x2-( M2-1 ) y2 ) 1 /2] 


+  g ( -Mt  -  M(M2-1  ) _1 [ -Mx  -  (x2-(M2-1 )y2) 1 /2 


2  1  /? 

for  x  >  0,  | y |  <  ( - 1  )  x. 


(16b) 


So  far,  the  analysis  has  been  carried  out  in  the  x,y,t  system. 
This  is  the  system  which  is  at  rest  in  the  surrounding  air.  The 
variables  x  and  y  appearing  in  the  last  equation  have  been 
defined  in  Eq.  (12);  they  have  been  introduced  to  simplify  the 
calculations.  Accordingly,  one  should  express  in  the  above 
formulae  x  and  y  by  x,  and  y.  To  obtain  the  corresponding 
expression  in  a  system  of  coordinates  in  which  the  source  is 
stationary,  one  must  transform  to  a  system  of  coordinates  which 
moves  with  the  source,  i.e.,  it  moves  with  the  dimensionless 
velocity  M  from  right  to  left.  Accordingly,  one  must  set 

x  -  x  +  Mt 

y  -  y 

4>  (  x  ,  y  ,  t )  =4>  (  x -Mt ,  y ,  t) 

But  x  and  y  are  identical  with  the  coordinates  x  and  y  introduced 
in  Eq .  (12).  The  right-hand  sides  of  Eqs.  (16)  have  therefore 
already  the  desired  form  if  the  step  from  x  to  x  is  not  carried 
out.  To  simplify  the  expression  further,  we  set 

Mg ( -ME )  »  g(E) 

One  will  remember  that  we  restricted  ourselves  to  one  meridian 

plane  with  coordinates  x  and  y.  In  general,  one  must  replace 

-2  -2  2  2  2 

ybyy  +  z  -  y  +  z  .  Then  one  has 

for  M  <  1 , 


g(t-(1 -M2)"1 [-Mx  *  (x2+(1 -M2) (y2+  z2))'/2]  ) 

r  2  /  ,  m2w  2  2  ,  1  /  2 

[x  +  (1-M  )(y  +  z  )] 


(  x  ,  y  ,  z  ,  t )  = 


(17) 


Cur 


and  for  M  >  1 , 

<t>  (  x  ,  y  ,  z  ,  t )  =  [x2-(M2-1  )  (y2  +  z2)  ]  _1  72 


g( t+(M2-1 )  1 [-Mx+(x2-(M2-1 ) (y2+z2) ) 1 /2] 


•g( t+(M2-1 )  1 [-Mx-(x2-(M2-1 )(y2+z2))1/2]) 


?  - 1  /? 

for  x  >  0,  y  <  (M  -1)  x,  otherwise  4>  =  0. 


The  different  appearance  of  the  fundamental  solutions  in 
subsonic  and  supersonic  flows  finds  an  explanation  in  Figures  la 
and  1b  which  show,  at  a  fixed  time  t,  perturbations  in  a  medium 
at  rest  caused  by  momentary  emissions  of  mass  from  a  point  moving 
from  right  to  left  with  the  dimensionless  velocity  M.  It  is 
assumed  that  at  t  =  0,  the  point  where  the  perturbation  occurs 
arrives  at  the  origin  of  the  x,y  system.  For  M  <  1,  the 
perturbation  generated  from  the  individual  points  lie  on  circles 
nested  within  each  other.  If  the  perturbations  started  at  x  -  * 
at  a  time  t  =  then  these  circles  cover  the  entire  x,y  plane. 

Only  one  circle  goes  through  a  given  point  x,y.  The  center  of 
this  particular  circle  gives  the  point  of  the  x  axis  from  which 
the  perturbation  started. 

If  M  >  1,  then  all  circles  lie  within  a  cone  through  the 
point  x  =  -Mt,y  =  0.  (One  will  remember  that  the  picture  is 
drawn  for  a  fixed  time  t.)  Inside  the  cone  one  finds  two 
circles  going  through  each  point.  More  familiar  is  the 
interpretation  of  these  figures  as  perturbation  introduced  at 
some  point  (x,y)  of  a  coordinate  system  in  which  the  air  moves 
with  the  dimensionless  velocity  M  from  left  to  right.  The 
perturbations  are  then  found  on  expanding  circles  whose  center 
moves  with  the  air  particles.  (Since  these  are  figures  referring 
to  just  one  instant  in  time,  both  interpretations  are  possible.) 
All  further  discussions  will  be  carried  out  in  the  x,y,t  system. 

For  an  oscillatory  source,  one  has 


^  V  *  /  V  V  V  V  V  V 


CSC! 


raws® 


Perturbations  at  a  Fixed  Time  Generated  by  a  Source 
Moving  from  Right  to  Left  with  Constant  Subsonic 
Velocity  (Fig.  la)  and  Constant  Supersonic  Velocity 
(Fig.  lb).  For  Supersonic  Flows,  the  Circles  Share  an 
Envelope  formed  by  a  Pair  of  Straight  Lines. 


g ( t )  -  exp(  ivt) 


Then ,  for  M  <  1 , 

^  r  2  /,  lJ(2w  2  2-.-1  /2 

4>  (  x  ,  y  ,  z  ,  t)  =  [x  +  ( 1  -M  )  (  y  +  z  ] 

(19) 

exp( iv(t-M(1  -M2)-1 [-Mx+(x2+(1 -M2)(y2  +  z2)  )1/2]) 
for  M  >  1 


4>  (  x  ,  y  ,  z  ,  t )  =  2[ x2+ (M2-1  )  ( y2  +  z2]~1/2 


(20) 


O  -1  P  -IP  P  PP1/P 

exp( iv( t-M(M  -1  )  x) )cos( v(M  -1  )  (  x  - (  M  - 1 ) ( y  + z  ) )  )  ; 

x  >  0,  | y |  <  (M2-1 )~1/2X. 


One  notices  the  occurrence  of  the  cosine  instead  of  an  exponential 
function  with  imaginary  argument.  The  difference  in  appearance  is 
the  reason  why  a  direct  transition  from  the  subsonic  fundamental 
solution  to  the  supersonic  fundamental  solution  is  not  possible. 

The  phenomenon  has  a  counterpart  in  the  general  unsteady  case. 

2  2  2  2 

As  (x  -  (M  -1)(y  +  z  ))  tends  to  zero,  (i.e.,  as  one  approaches  the 
Mach  cone  through  the  point  x  »  0,  y  »  0),  the  arguments  in  the  two 
expressions  in  Eq.  (18)  approach  each  other  and  one  can  develop  g 
with  respect  to 


One  obtains 


4>  (  x  ,  y  ,  z  ,  t )  =  2[x2-(  M2-1  )  (  y2  +  z2)]  1 /2[  g(  t-M(  M2-!  )  ]x 

♦  ( 1 /2(M2-1  ) _1 g'  '  ( t-M(M2-1) _1 x) (M2-1 ) “2(x2-(M2-1  ) (y2  +  z2)  ) 

+  ( 1 /24) (M2-1 ) _1 gIV( t-M(M2-1 )x) 

lyi 2  .  ,  -4  ,  2  lyi2  .  w  2  2,  .2 

(M  -1)  (x  -(M  -1 ) ( y  +z  ))  +.. 

This  situation  can  be  expressed  in  the  following  manner.  Let 

D  r  2  ,2  ,  w  2  2,  ,-1/2 

R  -  [x  - ( M  —  1 ) ( y  +z  )] 

At  the  Mach  cone  one  has  R  =  0.  Then  one  has,  from  Eq.  (18) 

4>(x,y,z,  t)=R_1  tg(  t+(M2-1  )  _1  (-Mx  +  R)  )+g(  t  +  ( M2- 1  )  _1  (-Mx-R)  )} 

For  fixed  t  and  x,  the  term  within  the  braces  is  an  even  function 
of  R. 

The  original  particular  solution  in  the  xy-system  is  valid 
throughout  the  xy-plane,  (which  is  fixed  with  respect  to  the 
undisturbed  air)  although,  because  of  the  introduction  of  the  6 
function  for  f,  the  individual  perturbations  are  restricted  to 
circles  and  large  regions  of  the  plane  remain  free  of 
perturbations.  Since  the  final  expression  in  the  xy-system  (in 
which  the  air  moves)  is  obtained  merely  by  a  coordinate 
transformation,  the  results  are  valid  in  the  entire  xy-plane,  in 
spite  of  the  singularity  which  occurs  at  the  surface  of  the  Mach 
cone.  This  means  that  the  partial  differential  equation  for  $  is 
satisfied,  even  as  one  passes  through  this  surface. 


SECTION  III 
UPWASH  CONDITIONS 


x  =  xL,  y=yL , z  = zL , t = tL/a 

i.e.,  x,  y,  z,  and  t  give  a  system  of  coordinate  corresponding  to 
x,y,z,t,  but  before  these  coordinates  have  been  made 
dimensionless.  Let 

A  A  A  A  A  A  A  A  A 

<j>(x,y,z,t)  =  aL<j>  ( x  ,  y  ,  z  ,  t )  =  aL<j>(x/L,  y/L  ,  a/L  ,  at/L  ) 

Let  the  deformation  and/or  motion  of  the  wing  surface  be 
given  by 

D( x  ,  y ,  t) . 

The  upwash  condition  is  then  expressed  by 


--  U 


3D  3D 

ax  at 


This  is  now  written  in  the  nond  imensional  form.  Let 


D(x,y,t)  =  LD( x/L , y/L , at/L) 


Then 


n  3D  -3D 

D  ,  =  a-r- 

X  3t  3t 


Since  U/a  =  M,  one  obtains 


>  =MD  +  D. 
z  x  t 


>.V.V.V.V.V.V.'AV  V  ■ 
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SECTION  IV 


THE  INTEGRAL  EQUATION  FOR  T IKE- DEPENDENT  SUPERSONIC  FLOWS 

General  time  dependent  supersonic  flows  over  a  wing  whose 
planform  lies  in  the  xy-plane  can  be  represented  by  doublets 
(oriented  in  the  z-direction)  distributed  over  the  planform  and, 
if  one  has  a  subsonic  trailing  edge,  over  part  of  the  wake. 
Fundamental  solutions  for  such  doublets  are  obtained  from  those 
for  sources  by  a  differentiation  with  respect  to  z.  (In  the 
present  context  a  normalization  of  the  sources  or  doublets  is  not 
needed.)  The  evaluation  of  the  upwash  due  to  this  (unknown) 
doublet  distribution  requires  a  further  derivative  with  respect 
to  z.  The  desired  integral  equation  for  the  doublet  distribution 
is  obtained  by  equating  the  upwash  so  found,  with  the  upwash 
prescribed  by  the  upwash  condition  Eq .  (21). 

Let  £)n  be  a  system  of  Cartesian  coordinates  in  the  plane  of 
the  planform;  £  corresponds  to  x  and  n  to  y.  Let  the  doublet 
strength  at  a  point  of  the  planform  (£,n)  be  given  by  h(t,5,n). 
Let  furthermore 


6  -  (M2  -  1)1/2  (22) 

To  obtain  the  potential  due  to  a  source  at  the  point  £  =  x,  n=y, 
one  must  replace  on  the  right  hand  side  of  Eq.  (18)  x  by  (x-£) 
and  y  by  (n~y). 

We  introduce 

R  =  [(x-£)2  -  B2z2  -  82(n-y)2]1/2  (23) 


and 


ret  =  6  2[M(x-£)  R]  (24) 

1  ,2 
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here  the  upper  and  lower  sign  pertain  to  the  former  indices  i=1 
and  i»2. 


The  symbol  ret  stands  for  retardation.  Roughly  speaking, 
it  gives  the  time  required  for  a  perturbation  to  propagate  from  a 
point  £ , n  to  a  point  x,y,z. 

In  the  coming  discussions  x ,  y  and  t  are  kept  constant.  In 
some  of  the  functions  that  will  occur,  these  quantities  are  not 
listed  as  arguments. 

The  potential  due  to  a  source  distribution  in  supersonic 
flow  is  then  given  by 

<{/s^(t,x,y,z)  =  j|(1/R)[h(t-ret1(x,y,C.n),^,n) 

+  h(t-ret(x,y,e,n),&,n)]d£dn  (25) 

2 

The  region  of  integration  comprises  all  points  of  the 
planform  and  sometimes  of  the  wake,  that  lie  within  the  forecone 
of  the  point  x,y,z,  for  only  then  is  R  real  and  x-£  >  0.  The 
boundary  of  the  region  consists  of  the  hyperbola  R=0, 

Eq.  (23)  and  of  the  leading  edge.  The  vertex  of  the  hyperbola 
lies  at 


(x-O  =  gz,  n  =  y  (26) 

for  z  >  0 

To  bring  the  portion  of  the  boundary  formed  by  this 
hyperbola  into  a  simple  form  we  introduce 

n  =  y  +  g  1 qn  ( 27  ) 

o  P  P  1 /? 

with  q(x-E,z,g)  -  [(x-£)  - B  z  ]  (27a) 

One  notes 

Then 


ia 

3  z 


-g2z 


d  “2.1/2 

R  -  q( 1 -n  ) 


(28) 
(29  ) 


The  hyperbola  R  =  0  is  then  transformed  into  three  straight 
lines  in  the  £n  plane 

n  “  1 

and  q  =  0 

Only  the  line  q  -  0  depends  upon  z.  The  condition  q=0 
becomes 

x-C  =  8  z 

Let  be  the  smallest  value  of  t  which  occurs  within  the  region 
of  integration.  Eq.  (25)  can  be  written  in  the  form 

x-8z 

<J> 3  ( t ,  x  ,  y )  =  |  (j^tjZ)  +  J2U,z))d£  (31) 

to 

where  after  substitution  of  the  retardation 

J1(2U,z)  - 

nupper  ^  ^ 

8_1|  (1-n2)_1/2h(t-8'2[M(x-C)iq(1-n2)1/2],t,y+61qn)dn  (32) 

n lower ^  ^ 


For  simplicity  x,y,  and  t  are  not  shown  as  arguments  of  J1  and 
J 2 .  For  those  values  of  t  for  which  the  upper  and  lower  limits 
of  n  are  formed  by  the  hyperbola  R-0  these  limits  are  _+1  .  For 
values  of  E,  where  one  or  both  limits  for  n  are  formed  by  the 


leading  edge  the  limits  depend  upon  £  and  z.  At  t .hose  limits  the 
potential  and  therefore  also  the  function  h  is  zero. 

To  obtain  the  potential  of  a  doublet  distribution,  with  the 
doublets  oriented  in  the  z  direction,  one  must  differentiate  this 
expression  with  respect  to  z.  We  consider  first  a  region  C1 <S<x 
where  is  sufficiently  close  to  x  so  that  the  limits  of 
n  are  _+  1.  Let  4)  ^ 3  *  ^  ^  and  <J>^d’^  be  the  contributions  to  the 
potentials  and  <f>^d^  from  this  region.  The  point  of 

departure  is 

x-Bz 

4>(S,I)  -  |  (J,  U.z)  ♦  J2(g,z)  )d£  (33) 


In  differentiating  this  expression  with  respect  to  z  one  obtains 
one  term  denoted  by  due  to  the  differentiation  with 

respect  to  the  upper  limit  and  a  second  term,  denoted  by  <J>^dIb) 
due  to  the  differentiation  of  the  integrand. 

One  finds 


(dla) 


=  -B[  J1  U-,  ,z)  +  J2U1  ,z) 


£=x-3z 


According  to  Eq.  (27a)  one  has  q  =  0  for  £  =  x-gz;  the  function  h 
in  Eq.  (32)  is  no  longer  dependent  upon  n.  Now 

+  1 

f  -2 i -1 / 2 
I  [  1  -n  J  dn  =  ti 


Therefore 

4>(dIa)  -  -2irh(  t-MB'1  z,  x-gz,y)  (3*0 

The  dependence  of  the  functions  J1  and  in  Eq .  (33)  upon  z 
enters  through  the  function  q  (see  Eq .  32).  As  q  appears  in  two 


of  the  arguments  of  the  function  h  one  obtains  from  each  of  the 
functions  and  J2  two  terms  for  their  derivatives  with  respect 
to  z  (superscripts  b  and  c).  We  denote  derivatives  of  h  with 


respect  to  its  first  second  or  third  argument  respectively  by 
hO),  ^(2)^  and  ^(3)^  Then  with  Eq.  (28) 

OJ1(b)/9z)  +  (3J2b)/3z)  = 

+  1 

=  ze"1q'1[lh(1 )(t-B'2[M(x-C)+q(1-n2)1/2] ,^,y+0_1qn] 

-1 

-h(1 )(t-6_2[M(x-C)-q(1 -n2)1 72 ] ,6,y+B_1qn)}dn  (35) 

and 

( 3 J 1 (c)/3z  +  (3J2(c)/3z) 

+  1 

= -zq~ 1 Jn(1 -n2)"1/2{h(3)(t-B"2[M(x-^)+q(1 -n2)172] ,C,y+8_1qn) 
-1 

+  h(3)  ( t-(T2[M(x-0-q(  1  -n2) 1  /  2J  +  qn)  tdn  (36) 

At  the  upper  limit  of  the  integral  (31),  one  has  q-O  .  The  factor 
q-1  in  the  last  Eqs.  (35)  and  (36)  might  cause  the  integrals  in 
Eq.  (31)  to  diverge,  but  because  of  the  special  form  of  these 
equations  this  does  not  happen;  the  reasons  are  different  in  the 
two  cases.  In  Eq.  (35)  the  two  terms  of  the  integrand  agree  for 
q=0,  except  for  the  signs.  Accordingly,  the  integrand  is  0(q); 
observing  that  for  q-0, 

x-£*8z 

and 
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+  1 

f  -2J/2.-  /0 

j  ( 1  -n  )  dn  =  it /2 

-1 

one  finds 

lim  (OJ  (b)/3z)  +  OJ?(b)/3z)  )  =  -znB’3h(  1  ’ 1  }  ( t-6_1  Mz,x-6z,y  ) 
q  -»0 


Regarding  Eq.  (36)  one  observes  that 


+  1 

j  n( 1 -n2) _1 /2h( 3) ( t-B'2[M(x-0 +q(1 -n2) 1 /2] ,5 ,y  )dn  =  0  (37) 

-1 

because  of  the  antisymmetry  of  the  integrand  with  respect  to  n. 

(  3 ) 

Notice  that  only  the  last  argument  of  h  differs  in  Eq.  (36)  and 
(37).  Subtracting  the  two  forms  of  Eq .  (37)  from  Eq .  (36) 
developing  the  expressions  with 
observing  that 

+  1 
f 
J 

-1 


one  obtains 


respect  to  the  last  argument  and 


-2  . - 
n  dn 


<i-;2)'/2  2 


lim  (  OJ  (c)/9z)  +  (3J  (c)/9z)  }  =  -ztt6  1  h (  3  ’  3  }  ( t -6  _1  Mz  ,  x -8  z  ,  y ) 

q  +  0 


Again  this  expression  is  bounded.  This  procedure  requires  that 
one  carry  out  the  integration  with  respect  to  n  first.  One  can 
now  carry  out  the  limiting  process  z  +  0.  Since  the  z  derivatives 
of  J1  and  contain  a  factor  z,  their  contributions  vanish  in 
this  limit.  One  therefore  obtains 


lim  =  lim  =  -2irh(t,x,y)  (38) 

z-"0  z*0 


3 

a 


m 
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Accordingly,  the  function  h  represents,  except  for  the  factor 
- 2 Tt  ,  the  potential  at  the  upper  side  of  the  wing. 


( d  I  ) 


In  order  to  obtain  the  upwash  due  to  the  expression  <{> 
one  must  form  the  derivative  with  respect  to  z  and  subsequently 
perform  the  limiting  process  z->0.  One  obtains  from  Eq.  (3*0 


lim  <j>(dIa)  =  2TT(MB"1h(1  )(t,x,y)  +  Bh ( 2 } ( t , x , y ) ) 


q  +  0 


or  in  more  conventional  notation 


lim  =  2tt  ( MB  'hh(t,x,y)  +  Bh  (t,x,y)) 


-1 


z  +  0 


Eqs.  (3b)  and  (36)  have  the  form  zF(z).  Then 


lim  •—  ( zF ( z  )  )  =  F ( 0  ) 
z-*0 


Taking  into  account  that  q|  -  x  -  £  one  obtains 

z  =  0 


x  + 1 

lim  4)(dIb)  =  [^47(  (|6‘1[h(1)(t-B"2[M(x-f;)^(x-E)(1'n?)1/2],C,y  +  B"1(x-^)n) 


h(1)(t-B-2[M(x-E)-(x-f;)(i-n^)l/^],f,,y  +  B”'(x-E)n)] 


-2,1/2 


■n(l  -n2)'1/2[h(3)(t-B’2[M(x-^)-(x-E)d  -n2)1/2]  (x-E)n: 


+  h^(t-8  2 [ M( x -5 ) - ( x -£ ) ( 1 _n2 ) 1 ]  ,  K  •  y  +  B  1 ( x -5 ) n ) ] } dn 


(  3 1 1 )  (  s  ) 

In  the  contribution  of  <J»  <}>  one  has  q/0.  The  upper 

limit  of  ,  namely  ^ ,  is  fixed.  In  differentiating  J1  and  J2 
with  respect  to  z  one  obtains,  in  principle,  contributions  from 
the  z -de r i vat i ves  of  the  limits  £lower  and  ^Upper’  While  these 
limits  are  fixed  in  the  F,n-system  they  will  depend  upon  z  in  the 
^n-3ystem.  These  limits  are  given  by  the  leading  edges  and  there 


h-0.  At  subsonic  leading  edges  a  square  root  singularity 
appears.  There  might  be  a  question  whether  or  not  this  has  an 
effect  on  the  second  derivative  with  respect  to  z.  This  has  been 
discussed  in  considerable  detail  for  the  steady  case  (Ref.  4) 
with  the  result  that  in  the  limit  z->0  there  are  no  contributions 
due  to  these  limits.  One  therefore  finds 

,  .  . ( sll )  - 

1 1  m  0  =0 

z-0 

Moreover,  lim  will  have  the  same  integrand  as  <j>fdI^;  only 

0z  z 

the  region  of  integration  is  different  Returning  to  the  original 
coordinate  n  one  obtains  the  following  expression. 

Let 


R  -  [(x-O2  -  82(n-y)2]1/2  (39) 

ret1  2  -  B‘2[M(x-U  R] 

Then 

4>d(  t ,  x  ,  y ,  z  =  0  )  =  2tt(M6  1ht(t,x,y)  +  Bhx(t,x,y)) 

*jj { h^ 1 ^ ( t-ret1 , £ , n ) -h^ 1 ^ ( t-ret2 , £ , n) 

-B2((n-y)/R)[h^\t-ret1,C,n)+h^^(t-ret2,C.h)]  }  ( 40 ) 

(x-t) 


The  region  of  integration  has  as  boundaries  the  leading  edge  and 
the  two  Mach  waves  through  the  point  (x,y)  into  which  the 
original  hyperbola  degenerates  for  z*0 .  Equating  this  expression 
with  the  upwash  found  in  Eq .  (21)  one  obtains  the  desired 
integral  equation  for  the  n(t,x,y).  From  Eq.  (38)  one  obtains 
the  potential  and  hence  all  quantities  of  physical  interest.  The 
integrands  in  Eq.  (40)  occur  with  a  retarded  time  argument.  In 
essence  they  are  known  from  the  results  for  preceding  time  steps. 


25 


y. 


If  data  up  to  a  certain  time  are  known,  then  Eq.  (40)  allows  one 
to  evaluate  h^(t,x,y).  This  is  the  essential  step  in  an 
integration  procedure. 


We  add  the  following  observation.  In  most  of  the  practical 
problems  the  upwash  computed  from  the  displacement  and  or 
deformation  of  the  wing  Eq.  (21)  can  be  written  in  the  form 


w(  t ,  x  ,  y )  =  l  c  (  t)  w  (  x  ,  y ) 

k  K  K 


where  the  number  of  functions  wk(x,y)  is  rather  small.  These 


functions  are  considered  as  known.  Whether  or  not  the  functions 


ck(t)  are  known  in  advance  depends  upon  the  nature  of  the 


problem.  In  aeroelastic  problems  the  determination  of  the 
aerodynamic  response  is  done  in  a  preparatory  phase  and  then  the 
deformation  are  not  yet  known.  In  the  aeroelastic  equations  the 


ck(t)  will  then  appear  as  dependent  variables.  It  is  practical 


to  express  the  functions  ck(t)  in  the  form 


C  k (  t )  =  c£0)  ♦  !*k(T)H(t-x)dT 


('ll ) 


Here  ck  denotes  the  derivative  of  ck  with  respect  to  its 
argument  and  H(t)  is  the  Hamilton  step  function 


H ( t )  =  1  ,  for  t  >  0 


H  ( t ) 


0  ,  f or  t  <  0 


One  observes  that  H(t-t)  =  1  in  Eq.  (41)  because  t>x, 
that  no  perturbations  are  present  at  t  •  0. 


We  assume 


Accordingly,  it  suffices  if  one  solves  the  integral 
equations  for 


w( t,x,y)  =  wk(x,y)  H( t-x) 
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separately  for  all  functions  w^(x,y)  that  are  of  interest.  Since 
t  does  not  occur  explicitly  in  the  formulation  of  the  problem,  it 
is  even  sufficient  if  one  solves  it  for 

w(  t ,  x  ,  y  )  =  wk  (  x  ,  y )  H  ( t )  . 

All  other  solutions  can  be  built  up  from  results  obtained 
with  these  expressions  for  the  upwash.  Since  for  t>0,  H(t)  is 
independent  of  t,  the  solution  for  h(t,x,y)  represents  the 
transition  to  a  steady  state  with  the  boundary  condition 
w  =  w  (x,y).  This  happens  within  finite  time. 

To  become  familiar  with  the  special  properties  of  the 
integral  equation,  we  consider  special  cases  in  the  next 
sections . 


SECTION  V 


SPECIAL  CASES 


A  nearly  trivial  example  is  obtained  for  a  two  dimensional 

steady  flow  in  the  xz-plane.  Then  w  =  w(x),  w  =  0,  h' 1 ^ =  h.  »  0 
(  o  )  or 

h  -  h  «  0.  The  integrals  in  Eq .  (40)  vanish  and  one  obtains 

w  =  w(x)  =  2Tt6hxt 
One  has,  according  to  Eq .  (38) 

<J>  =  -2tt  h(  x  ) 


By  definition  w  =  $  . 

Therefore,  from  Eq.  (40) 

<t> z  =  valid  for  z  =  0 

This  relation  is  a  solution  of  the  equation  for  two-dimensional 
steady  flow 


More  insight  is  obtained  by  the  one-dimensional  unsteady 
case.  Here  all  solutions  can  be  built  up  from  the  solution  of 
the  problem  with 


w  =  w(  t)  »  H(  t) 

From  the  partial  differential  equation,  which  in  this  case 
simplif ies  to 


zz 


tt 


0, 


one  expects  to  find 


z 


(42) 


We  therefore  set,  tentatively, 


ht  -  cH(t)  ,  hx  =  0,  hy  =  0 


where  the  constant  c  remains  to  be  determined.  The  integral 
equation  then  gives 


<t>z  =  H  (  t )  =  c  [  2n  Mg  H  (  t )  +  I ,  -  I2. 


where 


f  fH(t-g"2(M(x-U-R) 
^  (x-U2 


d^dn 


(44a) 


R  -  [(x-O2  -  e 2 ( n -y ) ^ ] l,d 

We  immediately  set  x=0  and  y*=0  (since  all  points  of  the  xy-plane 
are  equivalent)  and  consider  a  fixed  time  t.  The  region  of 
integration  lies  between  the  Mach  waves  through  the  point 

f>x  =  0  n  =  y  =  0 

i.e.,  between  the  straight  lines 

| n |  -  B'1^,  5<o 

As  a  further  boundary  one  has  the  curve  for  which  the  arguments 
of  the  functions  H  vanish 

t-g~2(-M5+R)  =  0  (45) 

The  more  general  case,  in  which  R  is  given  by 

D  r  ,  r .2  .2  2  2  ,  .  2i  1 /2 

R  -  l ( x-c )  -  B  z  -  8  (n-y)J 


2 1  1  /  2 


is  discussed  in  Appendix  B.  One  finds  for  both  signs  of  R 
in  Eq  .  ( ^a) 


U-(x-Mt))2  +  (n-y)2  -  t2  -  z2 
This  is  a  circle  with  center  at  E  =  x-Mt,  n=y  and  radius 

c . 2  2,1/2 

(  t  -z  ) 

In  the  present  case  z=0,  x=0,  y=0.  The  last  equation 
reduces  to 

2  2  2 

U  +  Mt)  +  n  =  t  . 

The  areas  of  integration  for  the  integrals  1^  and  are 
then  given  respectively  by  Figures  2a  and  2b.  Since  in  the  areas 
of  integration  H=1 ,  one  obtains 

T  T  f f d^dn 

:1  '  :2  =  ‘  )J-J2~ 

where  the  area  of  integration  is  the  circle 

U  +  Mt)2  +  n2  =  t2 


Setting 


5  +  Mt  - 

n  -  nt 


one  obtains 


Figure  2a.  Regions  of  Integration  for  z=Q  at  Different  Times. 

The  Boundaries  are  Formed  by  the  Mach  Lines  and  by 
Parts  of  Circles. 
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First  integrating  with  respect  to  n  one  obtains 


l1  1 2 


-2 


f 1  d^( 1 ~g2) 1 /2 


-1 


U-M)' 


The  integral  can  be  evaluated  in  terms  of  elementary  functions 
( Appendix  A ) . 


1,-12  -  2tt(1  -  (M/B  ) ) 

Substituting  this  into  Eq .  (44)  one  obtains 

H(t)  -  1  =  c  2tt 


Hence  c  =  (2tt) 

With  Eq.  (43)  one  finds 

ht  -  (  2-rr  )  H  ( t ) 


and  therefore 


4>  ^  *  -H  ( t ) 

This  is  indeed  the  result  anticipated  in  Eq.  (42) 

The  jump  of  at  t-0  encountered  here  is  typical  for 

p rob lems  whose  time  dependence  is  given  by  a  step  function.  This 

is  no  reason  to  abandon  such  an  approach;  it  is  the  key  step  to 

the  solution  for  general  time  dependence.  It  may,  however,  make 

special  measure  in  a  numerical  approach  necessary.  4>  and  its 

derivatives  are  functions  of  t,  x  and  y.  In  the  numerical  work 

one  will  therefore  generate  tables  which,  for  individual  points 

(x,y)  of  the  planform,  give  <J>  (and  probably  also  <j>t)  as  functions 

of  t.  <j>  and  <p  appear  in  the  integral  equation  with  the  first 
t  rj 


ffgy  .i.i  .*■  ■.-  •w±"y.,yr  ,r.u.,>v:-u.  -v-ft  -■-  v  .y.y.  -va -v./..* ^ 


«> 


i 

ft 
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argument  t-ret,  where  the  retardation  depends  on  x-e;  and  n~y.  In 
general  an  interpolation  will  be  necessary,  to  find  the 
quantities  in  question  for  these  arguments. 

As  far  as  the  determination  of  the  integrands  at  given 
points  ( E; ,  n  )  is  concerned,  this  does  not  generate  a  difficulty. 
One  will  always  have  t-0  as  first  value  of  the  independent 
variable  in  the  above  mentioned  tables  for  the  selected  points 
(x,y).  If  the  argument  t-ret  <  0,  one  always  obtains  zero;  for 
t-ret  >  0,  one  will  use  interpolation. 

In  this  example,  we  have  seen  that  the  vanishing  of  4>t  for 
t-ret-  <  01  defines  a  boundary  for  the  region  of  integration. 
Describing  the  situation  in  a  different  manner,  we  consider  the 
integrand  versus  n  at  a  fixed  value  of  E, ,  Figure  3.  At  t-ret  =  0 
the  integrand  jumps  from  zero  to  some  finite  value.  The  jump 
will  in  general  not  occur  at  one  of  the  chosen  points  (x,y).  If 
one  uses  interpolation  then  one  replaces  the  jump  by  a  ramp.  The 
error  so  introduced  will  be  small,  if  the  points  ( 5 .  n )  are 
closely  spaced.  It  may,  however,  be  preferable,  to  identify  in 
advance  the  point  where  the  jump  occurs  (here  the  value  of  n )  and 
modify  the  integration  formula  accordingly. 

Some  numerical  experiments  will  be  needed  to  arrive  at  a 
satisfactory  compromise  between  complexity  of  the  programs, 
computational  labor,  and  accuracy. 

In  the  next  section  we  treat  by  analytical  means  the  two- 
dimensional  problem  with  a  supersonic  leading  edge  and 
w(x,y)  *  H ( t ) .  In  this  case  one  can  represent  the  solution  by  a 
distribution  of  known  sources  and  find  the  potential  by  a  direct 
integration.  The  result  can  serve  to  test  the  accuracy  of 
numerical  approaches  where  the  same  problem  is  solved  by  means  of 
the  integral  equation  (which  uses  the  concept  of  a  doublet 
distribution) . 
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Actual  Integrand 


Boundary  Point  of  the  Region 
of  Integration  {t-ret)=0 


Figure  3.  Integrand  versus  n  at  a  Fixed  Value  of  The  Circles 

are  Points  at  which  the  Integrand  is  Tabulated  if 
the  Problem  is  Discretized.  If  the  Boundary  Point  of 
the  Region  of  Integration  is  not  Identified  then  the 
Jump  at  the  Boundary  will  be  Smoothed  out  over  one 


SECTION  VI 


THE  TWO-DIMENSIONAL  PROBLEM  WITH  A  STRAIGHT  SUPERSONIC 

LEADING  EDGE 

The  problem  can  always  be  reduced  to  one  in  which  the 
leading  edge  is  normal  to  the  flow  direction,  by  considering  only 
the  velocity  component  normal  to  the  leading  edge. 

If  the  potential  is  represented  by  sources  with  strength 
h(t,L,n),  then  according  to  Eq.  (38),  taking  into  account  that 

s  d 

4>z  (  t ,  x  ,  y  )  =  <(>  ( t ,  x  ,  y  )  =  -2irh(t,x,y) 

We  choose 

w(t,x,y)  =  2irH  ( t ) 

then 

h( t , 5 , n )  =  -H(t) 

and  from  Eq .  (25) 

fs)  [  [  —  p  o 

♦  (t,x)  =  -j]l/R[H(t-6  [M(x-L)+R]  )+H(t-6'"[M(x-L)-R]dE,dn 

R  =  [(x-L)2  -  B2n2]1/2 


The  areas  of  integration  can  be  reduced  to  those  points  of  the 
plan  form  for  which  the  arguments  of  the  unit  step  function  H  are 
positive.  Accordingly  one  writes 


0(3)(t,x) 


f  f  d£dn  f  [  dLrin 

JJ  ~R  +  JJ  R 


where  R1  and  R^  are  the  respective  areas.  The  regions  R1  an!  R, 
are  bounded  by  the  Mach  waves  through  the  point  n-0,  i  .  *■> .  , 


and  portions  of  the  circle 

(x-C  -  Mt)2  +  n2  =  t2  (48 ) 

This  circle  is  tangent  to  the  Mach  waves.  Tangency  occurs  at 
o 

- 1 

x-£  =  — j^—  t.  These  points  of  tangency  divide  the  circle  into 
two  parts;  they  constitute  portions  of  the  boundaries  of  R ^ ,  and 
R0.  Furthermore  the  regions  and  R?  are  restricted  to  the 
planf  orm. 

We  call  temporarily  basic  regions  of  integration  the  regions 

with  the  boundaries  given  by  Eq .  (47)  and  (48).  Their  shape  is 

seen  in  Figure  4,  case  a.  The  actual  regions  of  integration  R1 

and  R^  differ  from  these  basic  regions  on  account  of  the  fact, 

that  upstream  of  the  leading  edge  there  is  no  upwash  and 

therefore  no  sources  are  encountered.  For  fixed  t  the  shape  and 

size  of  the  basic  regions  of  integration  is  always  the  same.  If 

x  is  varied  then  the  entire  figures  shift  in  the  x  direction. 

For  different  values  of  x  different  parts  of  the  basic  regions  of 

integration  are  cut  off  by  the  leading  edge.  Thus  one  obtains 

different  cases  depending  upon  the  value  of  t/x.  They  are  shown 

in  Figures  4.  In  determining  which  configuration  applies  in  a 

specific  case  one  remembers  that  the  most  upstream  point  of  the 

circle  lies  at  £»x-(M+1)t,  the  point  of  tangency  between  the  Mach 

- 1  2 

wave  and  the  circle  lies  at  £=x-M  (M  - 1 ) t ,  and  the  most 
downstream  point  of  the  circle  lies  at  £=x-(M-1)t2.  The  leading 
edge  lies  at  £=0.  One  thus  obtains  the  following  cases, 

Figures  4. 


4.  Boundaries  of  Integration 
for  Fixed  z  and  Different 
Values  of  x  in  a  Supersonic 
Flow  with  Straight  Leading 
Edge  Normal  to  the  Flow 
Direction . 

For  each  Value  of  x  one 
has  two  Different  Areas 
of  Integrations.  The 
Shape  of  the  Areas  of 
Integration  is  the  same 
for  all  Values  of  x, 
except  that  Parts 
Upstream  of  the  Leading 
Edge  must  be  Excluded. 


Case  a 


t/x  <  ( M+ 1  ) 


Case  b  (M<-1)_1  <  t/x  <  M( M2 -1 ) ~ 1 


Case  c  M(M2-1)  1  <  t/x  <  (M-1)'1 


Case  d  ( M-1 )  <  t/x 


We  introduce 


fl  ft  -  1  ~ 

n  =  zr  ■ ,  then  n  =  B  n(x-£.) 

X  C,  f 

~2  1  17 

R  =  ( x-C ) ( 1 -n  ) 


The  integrals  in  Eq.  (46)  then  assume  the  form 


_  dCdn 

rrrr 

( 1  -n  ) 


The  integrations  with  respect  to  i  are  carried  out  first.  To 
de^rmine  the  limits  of  integration,  we  determine  the  points  of 
intersection  of  the  straight  lines  n=const  in  the  ^n-plane  with 
the  circle  that  forms  part  of  the  boundary  of  the  region  of 
integration.  Substituting  the  definition  of  n ,  Eq.  (49),  into 
the  equation  of  the  circle,  Eq .  (48),  one  obtains 

((x-S)-Mt)2  ♦  (x-£ )2n2/(M2-l )=t2 

This  equation  is  first  solved  for  (x-£),  but  the  result  is 
written  down  for  £ . 

Let  B,(;>  ■  WHilO-  ,  .:?).'(!  (50) 

M  -(1-n )  "  M  -( 1  -n‘~ ) 


Then  the  two  roots  for  £  are 


£ 1  =  x  -  (  B^  (  n )  +  B0(n)) 

£ 2  =  x  -  (  B1  (  n  )  -  B2  (  n ) ) 

These  are  the  lower  limits  for  £  ,  if  the  line  n=const  ends  on  the 
circle.  If  the  line  £;  =  const  ends  at  the  leading  edge,  £  =  0. 

Thus  the  lower  limit  is  given  by 

in  the  region  R^  by  max(0,£^(n)) 

and  in  the  region  by  max  ( 0  ,  (  n  ) ) 

The  upper  limit  is  always  £=x 

Then  one  obtains  along  a  line  n=const 

| d^  =  x  -  max(0,^1 (n) 

•=  min  (x  ,B1  (n)  +  B2(n))  in  the  region 
f 

j d^  =  min  ( x , B^ ( n ) -B^ ( n ) )  in  the  region 

Thus 


(S) (t,x)  =  -  8  ~ 1 [ 


n  j-  min(x,B1(n)  +  B,(n)) 

J - 


n  =  + 1 


dn  + 


n  =  -l 


2.1/2 
( 1  -n  ) 


n  =  -l 


min( x , B. ( n ) -B~ ( n  )  - 
- dn  ] 

"2,1/2 
( 1  -n  ) 


(51  ) 


The  transition  from  one  form  of  the  minima  to  the  other  occurs  at 
the  points  of  intersection  of  the  circle 


Case  b: 


( M  + 1 )  <  t/x  <  M  (  M  -  1  ) 


for  -  j  n 0  |  <  n  <  |  nQ  |  ,  min  (x.B^Bp)  =  x 


for  -1  <  n  <  -In, 


and  |  Dq  |  <  n  <  1  ,  min  (x.B^+B.,)  = 


for  -1  <  n  <  +  1,  min  (x.B^-B^)  =  -  B^ 


Case  c 


M(M2-1 )  1  <  t/x  <  (M-1 ) _1 


for  -1  <  n  <  1  min  (x.B^+B^)  =  x 


for  -1  <  n  <  - | nQ | 


and  | nQ |  <  n  <  1 


min  (x,B1-B2)  =  x 


for  -  |  nQ  |  <  n  <  |nQ|  min  (x,B1~B2)  -  B1  -  B2 


Case  d: 


(M-1)  <  t/x 


f or  -1  <  n  <  1 


min  ( x  ,  B^ =  x 


min  ( x , B^ -B2 )  =  x 

These  expressions  are  now  substituted  into  Eq.  (51),  the 
resulting  integrals  are  combined  in  such  a  manner  that  the  limits 

A  A  A  A  A  A 

are  either  n-_1  and  n=+1,  or  n  =  -n0  and  n=+n0.  One  obtains 
Case  a ,  t/x  <  ( M+ 1 )  1 

. <B, »B2»dn  *}  <B,-B2)dJ  , 

(  t ,  x  )  »  8  1  |  ~~2  i  /p  I  " p  1  /?  i 

_\  (1-n  T  2  ( i  -n2 ) 1  2 


.  *  .  '  »-»  »  .  v  .  ■  j-  w  v  t.  r.  <.  < 


<t>  ( t ,  x  )  =  -28 


-1  f  B1dn 

i  ""p  i  7p 

i1  ( 1 -n  ) 1  2 


Case  b,  (M+1)'1  <  t/x  <  M(M2-1)-1 


0  ( t, x  )  = 


(  B.  +B„ ) dn 


■1,  f  V£VB2 


J  , .  .2, 1 /2  +  J 
;  C i -n  )  i 


C  x-(B1+B2)dr, 


,,  -.2,1/2 

( i  -n  ) 


' _j  W)1/2' 


*(t,x)  --62 


-1  l,  I 


B1  dn 


J  77^27 T71  +  j 
_ ,  ( 1  -  n  )  - 


x-(B1 +B2)dn 

7r_’r2TT72’ 

( 1  -  n  ) 


Case  c,  M(M2-1)  1  <  t/x  <  (M-1)  1 


0 ( t,x ) --B 


•lif  xdn 


+1  D  1 0  * 

f  xd;  f  xd;  .  f  (BrB2)dn 

i  (1.fi2)f72an  J  M.fi2,1/2  J  ,,_ft2,1/ 


i  ( l ~n  ) 

-n 


.  ,  .2,1/ 

U  -n  ) 


0  ( t ,  x  ) 


Case  d, 


0  ( t ,  x  ) 


-B"1  1 2x  j - ^-f72 

^  (  1  -n  ) 1  72 


(M-1)  <  t/x 


+  1 

I  dn 

_2B  x  j  2  i  / p 

J  (1-n2)1  2 


f  (-x+B1-B2)dn 

i  — 


The  integrals  occurring  in  these  expressions  can  be  expressed  by 
elementary  transcendental  functions 

f  dn  .  ~ 

i - o~V7T<  =  arc  Slnfi 

Ml-n)  ' 


One  obtains  with  Eq .  (50) 


f  B2dn 


JU-nV/? 


=  t(M2-1)  [~d~-2  =  t(M2-1  )1 /2arctg  ^/2 


'  M  - 1  +  n 


( M  —  1  ) 


C  B2dn 


=  ^arctg  g 


f 


B,  dn 


J 


(i  -n  ) 


727TT2  ■  hzs 


dn 


( M  - 1  +  n  ) ( 1 -n  ) 


2,  f 72 


=  t8  arctg  (jj - %~T7o- 

6  (1-n2)1/2 


These  formulae  can  be  verified  by  differentiation.  Regarding  the  last 
expression,  it  is  shown  in  Appendix  B  how  it  can  be  derived  from  an 
integral  found  in  Ref.  3* 

One  then  obtains  the  following  specific  formulae 


for  t/x  <  (M+1 ) 


for 


for 


4>  (  t ,  x  ) 


<J>(x,t)  =  ~2tir  =  -2xit(t/x) 


(M+1 ) ~ 1  <  (t/x)  <  M(M2-1 )  1 


4>(x,t)  =  -2x[n(t/x)  +6  1  arc  sin  n 


( t/x )  arctg  ( |  - 


'0 


-2.1/2 

(i-n0) 


) 


( t/x) arctg  —g] 


M(M2-1 )  <  t/x  <  (M-1 )_1 


■2x[b  1  arc  sin  nn  +  (t/x)  arctg  ( 77 


6  n  ~  2,1/2' 
(  n0  5 


-  (t/x)  arctg  (-^)] 


mi 


for  (M-1)  <  t/x 


>.  \  •-  \  \1 ■. 


,  ^  ‘  .  '  .  *  -  '  . 


APPENDIX  A 


EVALUATION  OF  SOME  INTEGRALS 


The  integral 


!  ajLllTl 

l  U~  M)2 


2,1/2 


can  be  decomposed  into 


2M  M  -1  .  d £ 

*  U-M)2  o-?)1/2 


One  has 


f  d£ 

!  TTFP5 


5  »  arc  sin£. 


The  other  integrals  are  found  in  tables,  for  instance  Ref.  3, 
Formula  236,  3c 

f  dx  1  a2-ax  ,  ,  v 

J,  w  2  2,1  /2  =  ,  2  2,1  /2  rC  Sin  a  |  x  -a  |  *  1 01  I  >  3 

( x-a ) ( a  -x  )  (a  -a  )  1  1 

Theref  ore 


J  U-M)(  1  -f;2)172 


+  1 


1  1  -Mr 

T2aro  Sln  _M~r 
(  - 1  )  1  **  M 


6—  1 


I  r  .  I -K  .  1 +M ! 

— ? --T7p[arc  sinTrr-arc  sin  r-flJ 
(M  -1  ) 


o  1/2 
(M^-1  )  1 


/,  „•  .■  ^ 


Formula  236,  3a  of  Ref.  3  reads 


d£ 


-1 


I,  .  k ,  2  2  . 1  72 

J ( x-a )  (a  -x  ) 


( k -1 ) (a2-a2)  L  (x-a)k  1 


,  2  2.1/2  r 

[  -- — iTTT— (2k-3)ct  I - 


dx 


I,  .k-1  ,  2  2 

;  ( x-a )  (a  -x 


(k-2) 


dx 


I ,  .k-2,  2  2.1/2 

J (x-a)  (a  -x  ) 


Hence  for  k=2 


+  1 


f 


d£ 


u-m)2(i-s2)1/2 


1  -M 


M7T 


Fm2-1  )3/2‘ 


o-c2)1/2 


TT-mT 


d£ 


^  (C-M) ( 1 -C2) 1 /2 


Thus 


+  1 


f  d^(l -C2)1 /2 


^  (5-M)2 


v  (  -1 


2M 


,  [-1  * 


(M2-1  ) 


f/2 


(m2-i)1/2 


In  the  next  integral,  a  decomposition  is  carried  out,  so 
that  the  first  order  poles  are  displayed. 


( 


dn 


1  f 


[- 


^  (n2  +  62)  ( 1  -n2)  1  72  2lS  (n-iS  )  ( 1 -n2) 1  ^  ( n  +  i B  )  ( 1  -n2) 1  /2 


]  dn 


Theref  ore 


f _ dn 

J ( n2  +  B2 )  ( 1 -n") 


2Tf72  =  l  Im  J 


dn 


( n-iS ) ( l -n2) 1 72 


(a.  i : 


Now  one  finds  in  Ref.  3,  formula  236,  3c 

log 


(  _  _ dx_ 

( x-u ) ( a2-xc  ) 


-1 


d  2  2  w  2  2 .  1  1  / 

-ax  +  a  ♦  [  (  a  -a  ) ( a  -x  )  J _ 

x  -a 


J  ,  ~  2  2.1  /2  *  /  2  2.T72 

j  t  v  -r.  Wo  -v  1  (a  -a) 


2  .,2 


Specializing  for  a=1,  a  =  iB,  x  =  n,  one  obtains,  with  B  =*M  -1  , 


n-iB 


' ( n-iB  )  ( 1 

*  -r  log 

-  -R  log 

*  ~R[lOg[ 


~~2TU2  =  "M  x°8  1 
-n  ) 


P  1  /  P 

( -iSn+ 1 +M( 1 -n  )  ) ( n+iB  ) 


o  o 

n^  +  B  ' 


M  2,1/2  2,1  /  2 ,  u  2,1/2, 

Mn(M+(i-n  )  +iB(i~n  )  (M+(l-n  )  ) 

2  2 
n  +B 


n  +B 


The  last  term  in  the  bracket  is  a  constant  of  integration 
which  can  be  disregarded. 

Therefore,  from  Eq .  (A.1) 


APPENDIX  B 


REGIONS  OF  INTEGRATION  FOR  z/O 


Before  he  discussed  the  one-dimensional  time  dependent  flow 
for  z=0  (Section  IV),  the  author  investigated  this  problem  for 
z/0.  One  obtains,  as  expected,  the  results  usually  derived  from 
the  partial  differential  equation.  In  this  appendix  an 
intermediate  result  is  shown,  namely  how,  for  fixed  x,y,z,  the 
regions  of  integration  change  with  time.  In  Section  IV  the 
corresponding  regions  have  been  shown  for  z=0.  It  is  of  interest 
to  observe  the  relation  between  these  cases. 

It  is  assumed  that  the  time  dependence  of  the  upwash  is 
given  by  a  Hamilton  step-function.  The  arguments  are  given  by 


_  2 

arg1  )2  =  t  -  ret  =  t  -  B  (M(x-£)  _+  R) 


(  B  .  1  ) 


with 


R  =  [(x-S)2  -  B 2 ( n-y ) 2  -  82z2]1/2  >  0. 

The  areas  of  integration  at  fixed  x,y,z,t  are  characterized  by 
the  conditions,  that  R  is  real,  and  that  arg^  >  0,  arg.,  >  0.  In 
the  present  special  example,  all  points  of  the  £n~plane  are 
equivalent,  there  is  no  leading  edge.  In  the  general  case  the 
region  of  integration  has  the  same  boundaries,  except  that  points 
(£, n)  upstream  of  the  leading  edge  are  excluded. 

Part  of  the  boundary  of  the  region  of  integration  is  given 
in  the  £n-plane  by  the  hyperbola  R=0 

(x-O2  -  B2(n-y)2  -  B2z2  =  0  (B.2) 

In  the  interior  of  the  hyperbola  R>0.  Only  the  branch  for  which 
x~C  >  0  is  of  interest. 


51 


S3 


The  vertex  of  this  hyperbola  lies  at 


£  =  x-0z ;  z  >  0) 


n  =  y 


Its  asymptotes  are  given  by 


n-y  =  10  (x-O 


Other  portions  of  the  boundary  of  the  region  of  integration 


are  given  by  the  curve  along  which  arg .  or  arg„  vanish. 


Substituting  0  one  obtains  from  Eq.  (B.1) 


2  .,2,1/2 


(M  -i)t  -  M  ( x  -O  =  1  [(x-O  -  (M  -i)(n-y)  -  (M  -i)z  ] 


Squaring  both  sides  to  remove  the  square  root  and  dividing  by 


M  -1 ,  one  obtains 


M2t2  -  t2  -  2M(x-Ot  +  (x-O2  ♦  (n-y)2  +  z2  =  0 


( 5 - ( x -M t ) ) 2  +  (n-y)2  -  t2-z2 


(B.3) 


This  is  a  circle  with  center  at  Ox-Mt,  n  =  y  and  radius 

(t2-z  ^  .  One  surmises  that  the  hyperbola  R  =  0  is  tangent  to 


this  circle.  Elimination  of  (n-y),  from  Eqs.  (B.2  and  B.3)  gives 
the  values  of  £  (or  rather  (x-O)  for  the  points  of  intersection 
of  the  two  curves.  Actually,  this  process  eliminates  z  at  the 
same  time.  One  obtains 


(M(x-0  -  (M2-1  )t)2  =  0 


The  fact  that  one  obtains  a  double  root  for  x-£  indicates  that 
the  curves  given  by  Eqs.  (B.2  and  B.3)  are  tangent  to  each  other 


V 


■ 


.li.n  .1% 


(if  they  have  points  in  common  at  all).  The  points  of  tangency 
lie  at 


£ 


x 


The  corresponding  values  of  n  are  found  from  Eq.  (B.3) 


(n-y)2 


m2-i 

1  M2 


z 


2 


The  points  of  tangency  divide  the  circle  for  which  arg1  =>  0 
and  arg^  =  0  into  two  parts.  At  the  point  of  tangency  as  along 
the  entire  hyperbola,  R  =  0.  In  the  interior  of  the  hyperbola 
R  >  0.  If  one  travels  along  the  circle  starting  at  the  point  of 
tangency  in  the  direction  of  decreasing  (x-£),  then  arg^  (upper 
sign  of  R)  is  zero,  and  arg2  becomes  positive;  this  portion  of 
the  contour  is  therefore  the  boundary  of  the  region  there 
arg1  >  0.  The  portion  of  the  circle  for  which  x-E  is  greater 
than  the  value  for  the  point  of  tangency  is  the  boundary  of  the 
region  for  which  arg2  >  0. 


The  smallest  value  of  t,  for  which  points  of  tangency  exist, 
arises  if  the  circle  is  tangent  to  the  hyperbola  at  its  vertex. 
Then  (n-y)  =  o  and 


t 


M(M2-1  ) 


-1  /2 

z . 


This,  however,  is  not  the  smallest  possible  of  these  circles. 

2  2 

The  radius  is  given  by  (t  -z  ).  The  smallest  circle  therefore 
arises  for  t=z.  This  is  a  point  in  the  interior  of  the 
hyperbola.  There  R  >  0.  These  circles  in  their  entirety  are 
boundaries  of  regions  where  arg0  >  0.  One  thus  obtains  the 
sequence  of  boundaries  shown  in  Figure  5.  Those  boundaries  in 
Figure  5b  that  consist  of  full  circles  are  not  present  if  z=0. 
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Figure  5b.  Regions  of  Integration  for  z^O  at  Various  Times.  Fo 
Small  Time  the  Boundaries  are  Circles.  For  Greater 
Times  each  Boundary  consists  of  Part  of  a  Circle  and 
Part  of  a  Hyperbola. 
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APPENDIX  C 

DISCUSSIONS  RELATED  TO  FUNDAMENTAL  SOLUTIONS 

Fundamental  solutions  can,  of  course,  be  derived  directly 
from  the  partial  differential  equation  for  the  perturbation 
potential  in  a  system  of  coordinates  in  which  the  air  is  in 
motion.  Here  one  has 


—  (  M  —  1  )  +  0  +  0  -2M0  .  _  4> ...  =  0 

Yxx  Yyy  zz  Yxt  Ytt 


x  =  x,  y  =  y(M2-1)1/2,  z  = 


and  consider  perturbations  which  are  harmonic  in  time 


0 (x  ,  y , z)  exp  ( ivt) 


Then 


>  +  4> 

yy  zz 


2 i vM  %  v 

~2 —  ^x  +  ~2 — 
M  -1  M  - 1 


The  derivative  <(>"  is  suppressed  if  one  introduces 

?  7  f -i vM~  > 

<J>  =  $  exp  (-p — xj 

M  -1 


One  obtains 


>  +  <J>  +  0~~  -  — ^  0  =  0 

xx  vyy  (m2_, } 2 


Particular  solutions  are  obtained  as  follows.  Let 


D  2  “ 2  , “ 2  ~2, 
R  =  x  -  ( y  +  z  ) 


(C.  1  ) 


(  C  .  2  ) 


(C.3) 


(C.4) 


(  C .  5 ) 


(C.6) 


| 


V  V 

i  *  •  •  •  m 


and  assume  that  <p  depends  upon  R  only;  then 


This  is  solved  by 


<p  =  (1/R)  exp  (^~-)R  (C.7) 

M  -1 

Returning  to  the  original  coordinate  one  has 

p  p  P  P  P 

R  -  x  -  ( M  - 1 ) (y  *z  ) 

4>  =  (1/R)  exp(iut  -  --j---  _+  R) 

M  -1  M  -1 

These  are,  of  course,  the  particular  solutions  encountered  in  Eq. 
(20).  But  the  present  derivation  does  not  show,  which  linear 
combination  should  be  taken.  Besides,  it  is  not  clear,  whether 
the  continuation  of  these  particular  solutions  outside  of  the 
cone  R  >  0,  is  given  by  <>  =  0 .  We  shall  see,  that  this  is  not 
always  the  case. 

Eq.  (C.6)  has  no  physical  significance  because  of  the 

splitting  off  of  the  factor  exp(  x)  in  Eq.  (C.5).  Since  this 

M  -1 

factor  depends  upon  v  this  procedure  has  no  counterpart  in  the 
general  nonsteady  problem.  The  following  observations  may, 
nevertheless,  be  of  interest. 

Let 


v 


Then  one  has 


4> 


xx 


(  C .  8  ) 


58 


It  is  remarkable  that  in  this  equation  u  occurs  in  the  second 
power,  while  in  the  solution  Eq.  (C.Y)  it  appears  in  the  first 
power.  If,  however,  one  separates  in  Eq .  (C.7)  real  and 
imaginary  part  one  has 

j^rCOs(vR),  and  -^sin(vR)  (C.9) 


~2 

and  the  development  of  these  expressions  proceed  in  powers  of  v  . 

One  can  think  of  solving  Eq .  (C.8)  by  a  development  with  respect 
~  2 

to  v  .  To  the  lowest  order  one  then  has 


*RR  +  R  ^R"  0 


with  the  solutions 


<J»  =  ( 1  / R )  and  <j>  =  1 


<(>  =  ( 1  /R )  is  the  first  term  of  the  development  of  (1/R)  cos  vR,  <j>«1 
is  the  first  term  of  the  development  of  (1/R)  sin  vR  (after 
multiplication  by  v  1  ) . 

There  is  a  remarkable  difference  in  these  expressions  which 
will  be  explored  further  in  this  appendix.  The  solution  <j>  »  1/R 
is  real  only  within  the  Mach  cone.  But  it  can  be  continued 
outside  by  <J>=0,  and  the  expression  satisfies  the  differential 
equation  even  as  one  passes  through  the  Mach  cone. 

The  expression  <j>  =  1  holds  everywhere.  But  because  of  the 
nature  of  supersonic  flows  one  does  not  admit  perturbations 
outside  of  the  aftercone  of  the  point  xyz.  Therefore,  one  is 
inclined  to  set  <j>*0  outside  of  the  aftercone.  But  then  one  fails 
to  satisfy  the  flow  differential  equation  as  one  passes  through 
the  cone. 


We  discuss  the  case  v=0  further;  it  is  the  first  step  in 
development  with  respect  to  v.  The  function  $  is  the 


perturbation  potential, 
in  the  velocity  field. 


its  gradient  describes  the  perturbations 
Let  us  discuss  the  perturbations  in  the 
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N '  ■ 


mass  flow  vector  pw  (where  p  is  the  density) 
density  of  the  basic  flow 


Apw  =  pw  -  PQUex 


Let  pQ  be  the 


■+ 

w  =  Ue  +  Aw 
x 


Aw  =  <p  e  +  d>  e  +  d>  e„ 
yx  x  Yy  y  Tz  z 


One  has 


p  =  pQ  +  Ap 


P-  =  1/a; 
dp 


The  perturbation  in  the  pressure  Ap  is  found  from  Bernoulli's 
equat ion 


Ap  =  -poU<t>x 


Then 


Therefore 


Ap  =  -  P  0U4>X /a  =  -P0M$x/a 


pw  =  PQ(1  -M4>x/a)(ex(U^4>x)  +ey4>y +  ez4>  z ) 

^ -  =  e  <p  (1-M2)+e  p  +  e* 
pO  xYx  yTy  zTz 

The  potential  equation  expresses  conservation  of  mass.  If  it  is 
satisfied,  then  one  has 

d  i  v  (  p  w )  =  0 

and  for  the  linearized  form  of  the  potential  equation 


d i v ( A  pw ) 


=  0  . 

In  the  cylindrical  coordinates  x,y  this  equation  assumes  the  form 
g-|(yA(pw)x  +  gy(yA(pw)y)  =  0 

This  equation  is  satisfied  by  introducing  a  "stream  function"  Ay 
for  the  perturbed  flow 

y  A  (  p  w  )  x  =  2ttAi^ 

yA(pw)y  =  -2n  A»px 


Then 


At|)  =  2v  j  (  yA  (  pw)  xdy  -  yA(pw)dx) 

If  one  chooses  two  points  A  and  B  (Figure  6)  and  connects  them  by 
some  curve  and  evaluates  the  last  integral  along  this  curve,  then 
the  result  is  independent  of  the  choice  of  this  curve. 

Therefore,  Aip  is  a  function  of  x  and  y.  The  difference  A^g-AiJ;^ 
is  the  perturbation  in  the  mass  flow  passing  through  the 
axisymmetric  surface  generated  by  rotating  the  curve  AB  with  the 
meridian  plane  around  the  x  axis.  Since  Atp  depends  upon  x  and  y 
only,  one  can  draw  "stream  lines"  of  the  undisturbed  flow  (lines 
A^=const).  The  total  mass  flow  through  the  surface  AB  is,  of 
course,  the  mass  flow  due  to  the  unperturbed  flow  plus  the 
contribution  of  the  perturbation. 

This  is  applied  to  the  fundamental  solution  (Figure  7) 

4>  =  (x  -(M  -1)y^)  ,  x  >0 ,  |  y  |  <  (  M  - 1  )  x 


Then 


(M2-1 )y(x2-(M2-1 )y2)  3/2 


A(pw)  =  (M2-1)(e  x  +  e  y  )  ( x2-( M2-1 ) y2 )  3/2 

x  y 


This  shows  that  the  stream  lines  pertaining  to  the  perturbation 

are  radii  y/x  =  const,  although  the  pertinent  velocity  vectors 

have  a  different  direction.  For  simplicity  we  set  in  the 

2 

following  discussions  M  -1=1.  Let  us  determine  the  values  of 
as  functions  of  y/x  We  set  Ai^  =  0  for  y/x-0.  It  suffices  if  we 
carry  out  the  integration  along  a  line  x  =  const.  Let  y/x  =  a 


A\J/  ( a1  ) 


( x2-y2) 3/2 


_  ,  2  2.-1/2 
2irx(x  -y  ) 


=  2tt(  (1  -a2)'1  /2-1  ) 


As  a^-*1,  A  ip  ( a  ^  )■»«•;  the  mass  flow  due  to  the  perturbation  within 
the  Mach  cone  is  infinite. 

According  to  the  derivations  of  Section  I,  the  particular 
2  2  -1/2 

solution  <j>  =  (x  -y  )  is  valid  throughout  the  xy  plane,  for 

y/x  >  1  it  should  be  continued  by  <f>=0,  and  one  must  take  into 
account  the  infinite  jump  at  y/x=1. 

To  get  some  insight,  we  replace  the  jump  by  a  narrow 
transition  region  1 -e  <  y/x  <  1  (Figure  8).  There  we  set 


x[  1  -( 1 -e ) 2] 1 72 


1  -y/x 


(2-e  ) 


J _ r  1  _  _y  i 

1/2  372lx  2 ' 


(C. 10) 


At  y/x  =  1-e,  this  matches  the  original  expression 


,  2  2, -1/2 

( x  -y  ) 


At  y/x  -  1 ,  it  gives  0;  i.e.,  it  matches  the  outer  field  y  =  0.  We 
compute  the  perturbation  mass  flow  in  the  transition  region: 


»-•  vv-- JAvy  a' 
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y=x( 1 -c ) 


=  2Tre'3/2(2-e)"1/2[(  (1/2)(y/x)2-  (2/ 3)(y/x)3 


y/x=  1 


y/x= 1 -e 


One  obtains  for  e  small  as  mass  flow  in  the  transition  region 


-2,e-1/22-1/2 


The  perturbation  mass  flow  between  y/x  =  0  and  y/x=1-e  from  the 
original  expression  <j>  is 


2 tt  [  (  2e ) 


-1  /2 


The  total  mass  flow  between  y/x=0  and  y/x=l  is  therefore  -2n . 

The  physical  picture  for  the  perturbation  mass  flow  (not  for 
perturbation  velocity)  is  a  very  concentrated  inflow  in  the 
transition  region  and  an  outflow  for  y/x  <  1-e,  which  becomes 
very  large  as  y/x  approaches  (1-e).  The  total  inflow  exceeds  the 
outflow  by  -nr  (Figure  9). 


In  the  transition  region  the  potential  equation  (the 
equation  for  conservation  of  mass)  is  not  satisfied,  in  other 
words,  one  will  find  sources.  This  is  seen  in  detail  if  one 
computes  the  perturbation  mass  flow  through  a  conical  surface 
y=ax,  1-e  <  a  >  1 ,  extending  from  x^  to  Xg.  We  choose  x^  <  x^ 
then  the  following  expression  is  an  outflow  from  the  region 
y/x  >  a,  for  one  travels  around  the  region  in  the  positive  sense. 
The  potential  is  given  by  Eq.  (C.10).  One  has 


dy  =  adx 


V.VAjyA'AVV.VA.V.V, 


Figure  9.  Perturbation  Mass  Flow  Vector  in  Supersonic  Source  Flow 
at  Some  Radius  within  the  Mach  Cone.  The  Perturbation 
Mass  Flow  is  Finite  in  Radial  Direction;  the  Vector 
tends  strongly  to  Infinity  as  one  Approaches  the  Surface 
of  the  Mach  Cone,  The  Total  Mass  Flow  within  the  Mach 
Cone  is  Infinite.  Concentrated  Flow  Back  at  the  Surface 
of  the  Mach  Cone,  exceeding  the  Total  Outflow  by  zt. 
Therefore  the  Total  Configuration  Amounts  to  a  Source 
with  Negative  Strength. 


f  ■  !!■  L  »  V* 


H  ->  ;L ,« .-J  O . 


The  perturbation  mass  flow  through  this  section  for  M  -1 
then  given  by 


1  is 


X2  X2 
2ir  f  (_4>vdy  -  <J>ydx )  =  -2ua  j  (-a<j>x  +  <J>y)xdx 


J  '  Tx 

x. 


'1 


1 


2ira  e 


- 3 / 2 „ - 1  /2  f  t  ,  1  2a  s  1]  . 

2  j  (a  (— p  -  -p)  +  —  Jxdx 

1  X  X  X 

X1 


2ira  ( 1 -a  )  ( 1 +2a  )  e  ^22  ^^logCxp/x  ) 


For  a  close  to  one,  one  obtains 


6.2  1/2  e  3/2  it  log  ( Xp/x-j  )  ( 1  -a ) 


(C. 11  ) 


There  is  no  flow  through  the  line  a=1 .  In  the  transition  region 

-1  /2 

1 -a  <  e.  The  last  expression  is  therefore  log(x2/x1)  0(e  ). 

The  expression  (C.11)  is  the  outflow  from  the  region 


x1  <  x  <  x2 ,  a  <  y/x  <  1  ,  since  there  is  no  net  flow  through  the 


combined  cross  sections  x=x^  and  x  =  x2. 


The  expression  therefore 

represents  the  total  of  the  sources  within  the  region  just 
described . 


- 1 12. 

This  total  behaves  as  e  as  z+0, 


There  is  an  outflow  as  one  approaches  the  line  y/x-1-e  from 
above,  and  no  outflow  as  one  approaches  it  from  below.  At  the 
line  y/x=1-e  one  therefore  has  concentrated  sinks,  which  swallow 
the  total  of  the  sources  from  the  region  a  <  y/x)  <  1. 


The  sources  in  the  interior  of  the  region  and  the  sinks 
along  the  boundary  y/x=1-e  can  be  combined  to  form  doublets. 


Since  the  distance  between  the  sources  and  sinks  is  0(e)  the 

'/2>. 


total  doublet  strength  of  the  transition  region  is  0(e 
is  the  interpretat ion  for  the  source  solution. 


This 
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If  one  tries  to  use  a  similar  interpretation  for  a  potential 


<|>  =  1  ,  y/x  <  1 
0=0  y/x  >  1 

which  is  the  first  term  of  a  development  of  ^  sin  v  R,  one  finds 
that  the  flow  in  the  transition  region  through  cross  sections 
x  =  const  is  proportional  to  x9.  The  total  source  strength  in 
the  transition  region  is  obviously  not  zero,  these  solutions  are 
therefore  not  admissible.  At  least  for  flows  that  are  periodic 
in  time  this  shows,  that  only  the  first  of  the  expressions  shown 
in  Eq.  (C.9)  can  be  used  to  form  a  fundamental  solution. 
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SUMMARY 


The  report  derives  an  integral  equation  for  the  linearized 
supersonic  unsteady  potential  flow  over  a  wing.  Every  integral 
equation  formulation  for  a  problem  that  appears  originally  in  the 
form  of  a  partial  differential  equation  presupposes  the 
availability  of  a  fundamental  solution.  Such  a  fundamental 
solution  is  available  for  the  problem  at  hand  in  the  literature. 
It  is  rederived  here  to  show  its  particular  properties;  further 
discussions  are  found  in  Appendix  C.  The  integral  equation 
originally  obtained  requires  that  one  carry  out  a  limiting 
process  in  which  one  approaches  the  planform  from  above  or  below. 
This  formulation  is  brought  into  a  form  in  which  this  limiting 
process  no  longer  appears  and  one  works  solely  with  information 
available  at  the  planform.  Examples  which  can  be  treated 
analytically  bring  some  properties  which  have  a  bearing  on  a 
numerical  approach  into  sharper  focus. 
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